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Abstract 

One-dimensional model of a system where first-order phase transition occurs 
is examined in the present paper. It is shown that basic properties of the phe- 
nomenon, such as a well defined temperature of transition, are caused both by 
existence of a border between the phases and the fact that only in the vicinity 
of that border it is possible for molecules to change their phase. Not only the 
model is introduced and theoretical analysis of its properties is made but also 
the results of Monte Carlo simulations are presented together with the results of 
numerical calculation of the distribution of energy levels of the system. 

1 Introduction 

The prime purpose of the paper is to introduce a simple model of the system where 
first-order phase transition occurs. One of possible scenarios is presented here, which 
may lead to such phenomena as melting-freezing or evaporation-condensation. The 
task is mainly of didactic character. It does not discover any new types of systems in 
question but, basing on a very simple model, shows what are essential properties of 
the first-order transition. The model is meant to provide an introductory exercise for 
those students who plan to get involved in research of such phenomena, that is why 
numerous simplifications are applied in the model. 

Due to simplicity, our model allows for formal (mathematical) calculation of some ther- 
modynamic quantities (e.g. temperature of the transition, latent heat per molecule) 
as well as for the use of numerical simulation methods. Both approaches were taken 
in the present work. 

It is well known that one of the simplest attempts to classify phase transitions 
takes into account existence of latent heat. Hence, first-order phase transitions are 
those that involve latent heat (e.g. melting, evaporation) while second-order phase 
transitions do not (e.g. ferromagnetic to paramagnetic transition). 

There are a lot of models which describe phase transitions of both types. They are 
standard models such as Ising model, lattice gas model, XY model [l] but also some 
innovatory ones [2] . The initial point for most of them is Hamiltonian of the system 
for which partition function (or another statistical quantity) can be calculated. Then 
the partition function is examined and discontinuities or singularities, which may be 
responsible for the presence of the phase transition, are searched [5]. 

In our model the situation is opposite to some extent. We start constructing the 
model on the intuitive premise and only then the formal analysis may reveal that the 
partition function of the system has singularity at the thermodynamic limit . 

What features of crystal-liquid or liquid-gas transitions, apart from the latent heat, 
are most characteristic? It can be said that it is a well defined temperature (for fixed 
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pressure) of the transition. Below this temperature we have only phase I (e.g. crys- 
talline solid) while above there is only phase II (e.g. liquid). Both phases (generally) 
differ considerably from each other (symmetry, viscosity, compressibility, etc.), what 
means that a very little change of temperature of the system causes a dramatic change 
in the arrangement of molecules. Main purpose of our model is to explain intuitively 
how it is possible. 

A process of random walking at discrete time was adopted in our model in order to 
describe dynamics. Basic properties of melting are shown to be caused by the fact that 
there exists a border between the phases. Only near that border molecules may change 
their attachment to the phase. Far off the border they are blocked by their neighbours. 
Such a situation results in an exponential distribution of the energy spectrum of the 
system, the fact being a condition for the relation between energy and entropy to be 
linear so that phase transition may take place. 



2 Description of the Model 
2.1 The Basic Element of the Model 

The quantum system with a free hamiltonian Hq is a basic element of our model. The 
spectrum consists of n + m+ 1 (discrete, non-degenerated) levels, which are shown in 
Fig. 1. Such a system is called a molecule. 

The lower n+1 levels, numbered from to n, are called the levels of the first phase 
while the upper m levels (from n + 1 to n + m + 1) arc the levels of the second phase. 
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Figure 1 : Spectrum of a single molecule 

Energy of the k-th level of the molecule is given by: 

J Ai k for k < n 

" \ + A2(fc - n - 1) for A:> n, 

where A*2 = A12 + Ai n is energy of lowest level of phase II and A12 is energy gap 
between the phases, the existence of with is not necessary but added to make our 
consideration more general. What is more, our model is operational even if we assume 
Ai = A2 = A12. 

Neither the evident form of Hq nor its ,, origin" is needed. However, it is better for 
intuitive understanding to imagine that the states of Ho are the states of the motion 
of a molecule rather than the states of its internal electron excitations. 

In the following part of our work we examine behaviour of such molecules which, 
due to the interactions with external environment (thermostat) , change their state and 
jump at random to a neighbouring level (upper or lower). The chance of jumping has 
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been limited to transitions between the neighbouring levels only. Such a limitation 
is merely technical simplification and does not affect the appearance of the phase 
transition. 

2.2 Scheme of the System and Rules of Transition 

In the present paper we examine a system consisted of N molecules (as introduced 
in the previous subsection), which are numbered from 1 to and set into a one- 
dimensional lattice (Fig. 2). 
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Figure 2: Scheme of the lattice 

Additionally, it was assumed that individual molecules can interact with each other. 
The interaction does not change the structure of levels in Fig. 1 but causes that some 
of them become unavailable for the molecule, what means that there appear some 
restriction rules (which are essential for the phase transition to come into being): 

1. the i"* molecule may jump up from the top level of the first phase (i.e. n*'' level) 
onto the lowest level of the second phase (i.e. {n+ 1)*'* level) only if the left (i.e. 
{i — 1)*'*) neighbour of the molecule is currently in the second phase, 

2. the i*^ molecule may jump down off the lowest level of the second phase (i.e. 
{n + !)"■* level) to the top level of the first phase (n*'' level) only if the right 
{{i + 1)*'') neighbour of the molecule is currently in the first phase, 

3. exceptions: first molecule (i = 1) is always able to jump from n*^ level to (n+l)*'' 
level, last molecule (i = N) is always able to jump from (n + l)"* level to n*'' 
level, 

4. only that distribution of states was always chosen as the initial state of the 
simulation, for which all the molecules in phase II lie on the left to all of the 
molecules in phase I. 

The foregoing rules set up a border between the phases (analogous to a free surface) 
in the system. Transitions of molecules from the first phase to the second phase are 
possible only at the border of the phases. A molecule, which gathers enough energy 
to jump into a different phase, is not able to do it until its appropriate neighbour has 
done it first so the molecule is blocked by its neighbour and has to wait for its turn. 
It is shown below that such a property of the model guarantees that the first-order 
phase transition can exist. 

2.3 Limitations on Our Model 

• our model is one-dimensional while systems in nature are usually 3-dimensional, 

• molecules in phase II are distinguishable in our lattice while in gases and liquids 
they are undistinguishable, 
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• molecules in phase I have energy levels which are mutually independent in our 
model while in crystals their energy levels are common (collective fonon excita- 
tions), 

• distribution of levels in both phases of our model are uniform; there are regular 
gaps between the levels, 

• in our model the space, which is occupied by each phase is compact (i.e. it is in 
one piece) while in the real systems the phases may by mixed with each other 
as in fog, for example, when one phase is a suspension in another phase. Such 
an assumption is of a special kind because it is not connected with a character 
of local interaction between molecules but introduces some global limitations to 
possible configurations of the system. Although non-natural, the assumption is 
necessary for mathematical simplifications. Our model is only a simple demon- 
stration of process of growing (decline) a single drop of liquid during condensation 
(evaporation) . 

3 Theoretical Analysis of the Properties of the Model 

3.1 Statistical Properties of a Single Molecule 

Let us consider a sing le (free0) molecule defined in Fig. 1. A canonical partition 
function for such a system is given by: 

where An is contribution of the states, which belong to phase I: 

x-^ / ,x 1 — cxp(— /3Ai(7i + 1)) . , 

An = Eexp(-^A,fc) ^ ^ (2) 

while Bn^m is a contribution of states, which belong to phase II: 

n V^' o ( M\* ^(h ^\\ \\ exp(-/3A^2)(l - exp(-/?TOA2)) 
Bn,m^ 2^ exp(-/3(Ai2 + (fc-n-l)A2)) = — . . . (3) 

fc=n+l ^ ^^^^ P^^> 

In low temperature (i.e. for An > Bn,m) a free molecule can be found in the first 
phase more often than in the second. If to > n -I- 1 (i.e. number of levels in phase II 
is greater than the number of levels in phase I) the opposite situation is possible. In 
high temperatures (when A„ < Bn^m) the probability that a free molecule is in the 
second phase is greater than the probability of finding it in first phase. There is a 
temperature at which the levels of both phases can be occupied by molecules with the 
same probability: 

An{l3*) = Bn,m{l3*). (4) 

It will be shown that such a temperature is the one of the phase transition in our 
model (Fig. 2). This temperature will be marked by T* and its reverse by /3*. 

3.2 Statistical Properties of the Lattice 

Let us introduce Hamiltonian of the lattice (Fig. 2) in an evident form: 

H = (SfLiHo + H,nt (5) 



'^i.o. without limitations from section 2.2 
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where © is direct sum, Hq is a free Hamiltonian of single molecule, and Hint is an 
interaction Hamiltonian responsible for restriction rules as given in section 2.2. The 
way hamiltonian Hint acts can be imagined as that: it increases energies of ,, forbidden" 
eigenstates of hamiltonian ®^iHo by a certain very high value fl what results in an 
insignificant probability of forbidden eigenstates to be occupied. 
Partition function of the lattice is given by: 

N 

Zim=J2^n,mA:i-' (6) 

4=0 

where An and Bn^m are given by Eq. (2) and (3). The formula ([6]) can be written in 
a following form: 

N N 



i=0 1=0 

1 _ „JV+1 1 _ ( -l^N+l 

= e{(3~nA^^-^+e{(}*-(3)B:!,m ; I, (7) 

1 — a ' 1 — a ^ 

where a = Bn,m/ An, and is a Heaviside function. 

The function zi can be given in a form of two components, therefore a limit can be 
calculated for each of them as ^ oo. It can be seen that z; function gets divergent 
as N approaches infinity (due to and B^). However, with (3 tending to (3* there 
appears a further divergence (divergence of the sum of geometric series) . 
Averaged energy of the lattice per single molecule is given by: 

-<^--«(-^-^|-^^)^ (S, 

For — > oo wc obtain: 

E{p) - e{p - (3*) {El) + 6{p* - 13) {En), (9) 

where {Ei) and {En) are average energies of molecule in the first and the second phase 
respectively. 

For phase I (when levels of phase H are unavailable) averaged energy is given by: 

- - Ai ^^^-^ - —^^^-^-—^ J (10) 

while in phase H (when phase I is unavailable) it is given by: 

= = + A, ( I ^^IL^') . (11) 

of3 \exp(pA2) — 1 exp(pA2m) — 1/ 

It can be seen that for iV ^ cxd an energy leap takes place at temperature T ~ T* (see 
eq. [9|), what let us define latent heat per single molecule in the following form: 

CT= lim_E{P)~ \im^EiP) = {Eii)^p,^-{Ei)^p,y (12) 
In the state of equilibrium a probability of finding k molecules in phase II is given by: 

nfe AN~k 

Pk = " . (13) 
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It is worth noticing that: 

Pk+i=aPk, (14) 

what allows (when taken into account that probabilities are normalised to unity) for 
eq. ()13|) to be written in a following form: 

where a = i?„^™/A„. 

With the use of the formula ([TS]) an average number of molecules occupying phase II 
can be calculated in the form: 



N 



(1 + iVa^+i -a^(l + Ar)) 



Nn^^kP, = (16) 

and average number of molecules occupying phase I is defined &b Ni — N — Nn. 
An approximate formula for energy of the lattice per single molecule as a function of 
temperature is given by: 

E{I3) = [Nj {Ei) + Njj {Eii)) N~\ (17) 

where averaged energies are given by eq. (fTU)) and eq. (|lip . Formula ([TT]) is an estima- 
tion as it was not considered that near the border between the phases molecules had 
not enough time to reach the equilibrium state in their current phase (i.e. molecules 
did not "forget" yet that some time ago they had occupied the other phase). 

3.3 Connection with Random Walking 

It is convenient to use Markov chains to make it possible to understand why phase 
transition occurs in our model. The states of the process are numbered to iV and 
they define the number of molecules occupying phase II (or in other words a position 
of the phase border in our lattice). 

That process is a random walking (Fig. 3), with probabilities of jumping right and left 
given by: 

Bn.rn /l o\ 

" " TTb — ' ( 

where a - probability of transition of the next molecule to phase II, P - probability of 
transition of the next molecule to phase I. Probabilities (HH]), are chosen so that 
a + P = 1 and a stationary distribution of Markov chain would be in accordance with 
the equilibrium distribution p^ - P^ . 

Referring to well known properties of random walking |1] we get the following facts: 

• for T < T* (i.e. for a < [3) stationary distribution of the process is exponential 
(left plot in Fig. 4), what means that in the state of equilibrium most of the 
molecules occupy levels of phase I. Only at the left edge of the lattice there are 
some fluctuations leading to the fact that some molecules occupy levels of phase 
II. With T approaching T* fluctuations grow stronger and stationary distribution 
tends to become uniform, 

• for r > T* (i.e. for a > 0) the situation is opposite (right plot in Fig. 4). 
Majority of the molecules are in phase II except for fluctuations at the right 
edge of the lattice which cause that some molecules occupy levels of phase I, 
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Figure 3: Markov chain of transition of molecules between the phases 

• when T = T* , it is a case of symmetrical random walking (a = /?) , which cor- 
responds to a uniform distribution of probability p„ = Large fluctuations 
of energy are produced during the phase transition, what is apparent if we re- 
mind an example of the system that consists of the mixture of water and ice, at 
temperature 0°C, and set in a thermostat containing water and ice. A canonical 
distribution is used there, and fixed temperature in the zone of phase transition 
(with latent heat) does not establish the energy of the system. 

For a very large lattice (when N 00) the effect of edge fluctuations on all the 
properties of the system is negligible, even if T « T* (i.e. a « /?). With the increase 
of temperature of the system it can be noticed that nearly all the molecules change 
their phase simultaneously when temperature is exceeding T*. Such a behaviour could 
be spotted in practice, in numerical simulations described in the next sections. 

Pn Pn 
0.004. 0.004 

O.OO3I 0.003 

0.002 \ 0.002 

0.001 \ 0.001 

1000 2000 3000 4000 50oB 1000 2000 3000 4000 50oB 

Figure 4: Stationary distribution for asymmetrical random walking N = 5000 (left 
figure: a = 0.501, (3 = 0.499, right figure: a = 0.499, P = 0.501) 

In our model, a transition of a molecule from phase II to phase I is possible only 
if it joins itself to another molecule that is already occupying phase I (with except 
for the first molecule, which may always jump to phase I). It is in correspondence to 
real solidification, where molecules most often join an already existing crystal or begin 
to crystallise on microscopic grains of dust. In our model the right end of the lattice 
plays a role of a grain, because molecule i — N can always jump to phase I. Had that 
rule been eliminated from our model, all molecules being converted to phase II in high 
temperature, then after another fall of temperature molecules would not be able to 
return to phase I what would produce a supercooled phase II. 

3.4 Examples 
3.4.1 Example 1 

Let n = 1, m = 3, Ai = A2 = A12 = 1 then: 

• temperature of transition, calculated according to Eq.Q with the use of a nu- 
merical method (Fig. 5 - left plot) is T* = 6.08152, 
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Figure 5: Calculation of transition temperature T* . Left plot: example 1 (section 
3.4.1), right plot: example 2 (section 3.4.2) 




Figure 6: Amount of molecules occupying phase II {x — Nn/N) as function of tem- 
perature (curves for N = 20, 100, 1000). Left plot: example 1, right plot: example 2 

• average energy of the molecule in phase I for T = T* calculated according to 
Eq.dini) is {Ei)^^p,) = 0.458984, 

• average energy of the molecule in phase II for T — T* calculated according to 
Eq.mi]) is {Eii)(^p,) = 2.89087, 

• latent heat per single molecule calculated basing on Eg. p^ is ct = 2.431886. 
3.4.2 Example 2 

Let n = 2, TO = 20, Ai = A2 = A12 = 1 then: 

• temperature of transition, calculated according to Eq. ^ with the use of numer- 
ical method (Fig. 5 right plot) is T* = 4.3602, 

• average energy of the molecule in phase I for T = T* (calculated according to 
Eq.dlOl)) is (i;/)(/3.) = 0.848427, 

• average energy of the molecule in phase II for T = T* (calculated according to 
Eq.lHI])) is {Eii){T*) = 6.67352, 

• latent heat per single molecule calculated basing on Eg. p^ is ct — 5.825093. 

The plots in Fig. 6 depict the amount of molecules in phase II {x = Nu/N) against 
temperature. The left plot represents the example 1 while the right - example 2. The 
plots in Fig. 7 depict average energy of the system per single molecule (Eq. [8|) as 
function of temperature (right plot for example 1, left plot for example 2). We can see 
that with increasing N the curves become more and more steeper near the temperature 
of the phase transition. 
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Figure 7: Average energy of the lattice per single molecule as function of temperature 
(curves for N = 20, 100, 1000). Left plot: example 1, right plot: example 2 

4 Numerical Simulation 
4.1 Monte-Carlo Simulation 

In order to verify our consideration in practice we created a computer program (sim- 
ulation in Java language). We examined lattices, which were described in examples 1 
and 2 (last section) with the use of Monte Carlo algorithm[5]. In our algorithm, each 
molecule changes its position at random to higher or lower at each step and probability 
of such a leap depends on temperature as follows: 



A 
B 



exp(-/3A) 
1 + exp(-/3A) ' 
1 

1 + exp(-/3A) ' 



(20) 
(21) 



where A - probability of jumping up to a neighbouring higher level (if transition is 
permitted), B - probability of jumping down to a neighbouring lower level (if transition 
is permitted), A - energy gap between the levels. 
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Figure 8: Average energy of the lattice per molecule as function of temperature for 
N = 1000 (left plot: first example, right plot: second example) 

The results of our simulation are shown in figure 8. Left plot is drawn for the first 
example (section 3.4.1), right for the second one (section 3.4.2). The solid line is added 
for comparison - it is a theoretical curve ([5]), formerly shown in Fig. 7. Fat dots are 
end products of numerical simulation. Each of them represents (averaged over time) 
energy of the lattice (per a single molecule), for 500000 time steps. 
Simulation in the form of Java applet is available at www... (soon) 
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4.2 Numerical Computation of Distributions of Levels and of 
Other Thermodynamic Properties 

For small lattices (i.e. for small N, m and n) it is possible to calculate energy for 
each state of the lattice directly (numerically) . It was done with the use of a computer 
program for the system described in the first example (section 3.4.1) and for iV = 20 
molecules. We have obtained: 

- 81 energy levels of the system (first level: = 0, last: E = 80), 

- 10458 256 051 states of the system of aU, 

- degeneracy c?„ of energy levels (i.e. number of states of the system with the same 
energy) - see appendix A. In Fig. 9 logarithm c?„ is plotted against energy £'„ 




Figure 9: Energy levels degeneracy dn (in logarithmic scale) as function of energy 

In the Fig. 9 a linear fragment of the curve can be noticed, which corresponds to 
exponential increase of level degeneracy. If degeneracy d„ was identified with density 
of energy levels p{E) and S = k In p{E) then S{E) would be linear and responsible for 
phase transition occurrence. 

It can be proved intuitively in a following way that density of levels increases expo- 
nentially: 

• Let Di be amount of states of the system, which are available when border 
between the phases is located near i-th. molecule (i.e. Nu = i). Then Di+i = 
Di'm/{n + 1), what means that Di increases exponentially with respect to the 
number of molecules occupying phase II (iV//), 

• refering to Eq. (fTT]) we can see that energy of the system during the phase 
transition [T — T*) is linear function of A^//, 

• both above-mentioned facts joined togehter let us conclude that during the phase 
transition density of energy levels should increase exponentially with respect to 
energy of the system. 

When numbers dn are known, they can be used at calculation of probability distribu- 
tion of levels to be occupied at fixed temperature with the use of the formula: 

Pn(/3) = ^exp(-/3i?„). (22) 

In Fig. 10 the distribution is plotted for several temperatures, the temperature 
of the transition being also included. It is worth noticing that the distribution for 
T = T* gets flat (as it was expected) since there are symmetrical random walking and 
large energy fluctuations then. 
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Figure 10: Distribution of probability of level occupation for temperatures T=2, 4, 5, 
T*, 7, 9 



To compare numerical results with theoretical prediction (according to section 3) 
an average energy of the system (per a single molecule) was calculated as follows: 



(23) 



and is shown versus temperature (fat curve in Fig. 11), where Pn are given by Eg. ((22)) 
and energies are equal to En — nA. The hardly visible thin curve in Fig. 11 (which 
almost merges with the thick one) represents the theoretically computed energy (as in 
Eq. (IH)) and was alredy presented earlier in the left plot of Fig. 7. We can see that 
numerical results and theoretical expectation of the last section are in accordance. 




Figure 11: Average energy of the system per single molecule versus temperature (com- 
parison) 
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Appendix A. Degeneracies of energy levels of the sys- 
tem 

d„={l, 20, 191, 1160, 5037, 16683, 43969, 95191, 173988, 275785, 390079, 507601, 
626181, 750774, 889367, 1049453, 1237218, 1458361, 1719007, 2026235, 2388372, 2815232, 
3318382, 3911457, 4610529, 5434542, 6405826, 7550702, 8900195, 10490875, 12365848, 
14575924, 17180994, 20251653, 23871113, 28137458, 33166302, 39093922, 46080951, 
54316731, 64024444, 75467160, 88954965, 104853349, 123592989, 145680735, 171710383, 
202368411, 238422447, 280667788, 329792568, 386115323, 449167565, 517155934, 586442932, 
651295360, 704192361, 736884109, 742127143, 715675486, 657851392, 574032176, 473729446, 
368478879, 269241736, 184195595, 117577731, 69769278, 38326383, 19398465, 8996217, 
3797379, 1447096, 492785, 147992, 38513, 8475, 1521, 210, 20, 1} 

Numbers listing in the bracket are degeneracies of energy levels of the lattice (from 
E = 0to E = 80) - see section 4.2. 
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